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ствует уменьшению В. Кинетика фильтрования 
резко осложняется при переходе к режиму вы- 
соких селективностей, так как концентрация тем 
сильнее изменяется вдоль потока, чем выше $, 
осадок растет неоднородно, ], а следовательно, и 
АЮ/АЬ в каждый момент времени являются убы- 
вающими функциями 7 (расстояния от входа в 
канал). Это уточнение теории может стать пред- 
метом отдельного сообщения. 
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Введение. Применение мембран из коллоидных 
частиц или полиэлектролитов (называемых так- 
же динамическими мембранами) для обессоли- 
вания воды методом обратного осмоса [Л1, 2] 
повышает интерес к исследованиям особенностей 
формирования, структуры, кинетических свойств 
таких мембран. Динамическую мембрану (ДМ) 
получают фильтрованием коллоидного раствора 
мембранообразующего компонента (МОК) через 
мембрану с крупными порами (она называется 
подложкой динамической мембраны). Как пра- 
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вило, получаемые таким образом ДМ находятся 
в динамическом равновесии с потоком мембрано- 
образующего компонента и разрушаются при от- 
сутствии в фильтруемом растворе МОК (хотя 
известны мембраны, в которых МОК коагули- 
рует необратимо — такие ДМ, однажды сформи- 
рованные, в дальнейшем не разрушаются при 
снятии давления или фильтрации раствора, не 
содержащего МОК). 

В дальнейшем будем рассматривать только 
истинные динамические мембраны, в которых 
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коллоидные частицы не’образуют прочных коа- 
гулятов и удерживаются на конечных расстоя- 
ниях друг от друга только под действием потока 
жидкости, который прижимает частицы к под- 
ложке. Среднее расстояние между частицами 
(и объемная доля частиц) в такой мембране 
определяется из условия равенства гидростати- 
ческого давления в данном слое мембраны и 
расклинивающего давления коллоидных частиц. 
Зависимость расклинивающего давления от объ- 
емной доли частиц (уравнение состояния дис- 
персии) в этом случае определяется величиной 
сил отталкивания между частицами, в частности 
сил электростатического отталкивания. 

Таким образом, ДМ можно назвать концент- 
рированной дисперсией, стабилизированной от- 
талкиванием перекрывающихся двойных слоев. 
Теория строения и кинетических свойств таких 
систем должна учитывать взаимодействие со- 
седних частиц: гидродинамическое, электричес- 
кое, концентрационное взаимодействия перекры- 
тых двойных слоев, причем при разработке 
математического аппарата теории следует отда- 
вать предпочтение методам, отражающим в мак- 
симальной степени физические особенности рас- 
сматриваемой системы. 

Ячеечный метод, предложенный для учета 
гидродинамического взаимодействия частиц в 
дисперсии [3], может быть применен и для вы- 
числения уравнения состояния дисперсии, стаби- 
лизированной отталкиванием двойных слоев [4, 
5]. Результаты этих исследований будут приме- 
нены в данной работе для изучения гидродина- 
мической проницаемости ДМ. 

Постановка задачи. Скорость фильтрации 
жидкости ий через динамическую мембрану бу- 
дем вычислять, как функцию приложенного дав- 
ления Р при условии, что М — число частиц 
МОК над единицей площади подложки посто- 
яННО: 


Ч=0 (Р) езезааи (1) 

Физически выражение (1) означает, что при 
изменении давления над ДМ изменяются тол- 
щина, пористость, число частиц на единицу’ пло- 
щади и т. д. Задача о нахождении полностью 
стационарных свойств ДМ должна исходить по- 
этому из концентрации МОК в исходном раст- 
воре, физико-химических свойств МОК и свойств 
течения раствора над мембраной (которые опре- 
деляют условия отрыва частиц от ДМ и перехо- 
да их в раствор; без последнего процесса ника- 
кое стационарное состояние истинно-динамичес- 
кой мембраны невозможно). Поскольку эта пол- 
ная задача сложна, ограничимся изучением 
свойств ДМ в неполностью стационарном со- 
стоянии. Данное состояние характеризуется 
именно постоянством величины №. Известно [1, 
2], что время формирования стационарной ДМ 
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Рис. 1. Эскиз временных диаграмм скорости, давления и 
числа частиц в ДМ. 


составляет несколько часов. Таков же порядок 
времени, необходимого для разрушения ДМ. Та- 
ким образом, процесс наслаивания новых частиц 
на подложку или отделения частиц от ДМ про- 
текает сравнительно медленно; между тем, про- 
цесс установления скорости при изменении внеш- 
него давления — очень быстрый, время релак- 
сации здесь порядка времени распространения 
звука в жидкости поперек ДМ. 

Общий характер зависимости И от времени 
при скачкообразном повышении давления в мо- 
мент & изображен на рис. 1. При #<& скорость 
имеет полностью стационарное значение, соот- 
ветствующее давлению Ру. После повышения 
давления скорость устанавливается за период ла 
формирования гидродинамического режима (он 
имеет порядок микросекунд); при этом геомет- 
рия мембраны еще не успела измениться. Далее, 
за период т» происходит релаксация упругих 
свойств мембраны (этот период — порядок вре- 
мени распространения звуковой волны в твердо- 
образном теле ДМ по толщине мембраны); за тз 
происходит релаксация пластических свойств 
мембраны и за т. — наслаиваются новые части- 
цы на ДМ. Порядок времени т. — несколько ча- 
сов, о временах релаксации то и тз судить край- 
не трудно и соотношение этих времен, представ- 
ленное на рис. |, произвольно. Можно, однако, 
предположить, что скорость звука в теле мем- 
браны составляет не менее нескольких метров в 
секунду, т2 — миллисекунды. Если ДМ состоит 
из отталкивающихся частиц, то пластических 
деформаций в ней не возникает (тз= оо). Таким 
образом, в течение довольно большого проме- 
жутка времени т.—1. мембрана сохраняет по- 
стоянное значение М, и если к мембране при- 
кладывать импульсы давления длительностью 
намного больше т, но намного меньше та, то из- 
менение числа № несущественно для вычисления 
скорости, и условия отвода частиц от ДМ несу- 
щественны. 
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Рис. 2. Эскизы профилей объемной доли (1) и давления (2) 
в ДМ при низких (а) и высоких (6) давлениях. 


Следовательно, гидродинамическая характери- 
стика (1) отражает реакцию уже сформировав- 
шейся при определенном давлении ДМ (именно 
этому давлению соответствует заданное М№) на 
кратковременные изменения давления. 

Метод решения. Поставленная задача решает- 
ся посредством конкретизации формул из [2] 
применительно к конкретному уравнению состо- 
яния дисперсии, полученному в [5]. Указанное 
уравнение состояния отражает зависимость рас- 
клинивающего давления в дисперсии от объем- 
ной доли частиц в предположении, что поверх- 
ностный потенциал частиц мал, а толщина двой- 
ного слоя сравнительно велика. 

Квазистационарное состояние ДМ (при Ё>Ь) 
характеризуется равенством 


Р(х) = Рак (@(х)), (2) 


где Рав= Рав (9%) — расклинивающее давление в 
дисперсии; Р(х) — давление в ДМ на расстоя- 
нии х от подложки (рис. 2); «(х) — объемная 
доля частиц в ДМ на расстоянии х от подложки. 
Полное число частиц над единицей площади ДМ 
равно 


1 Е р | |2 
х 
= | п(х = \ = -=\ аха (х), 3 
м- [29д- [5р-в| 428. © 
0 0 
где п — плотность частиц; 6 — радиус ячейки; 
а — радиус частицы; Ё — толщина ДМ. 
Скорость фильтрации жидкости через ДМ 
определяется уравнением 


= — Га 4 = соп5, (4) 


где Гл. — гидродинамическая проницаемость диспер- 


сии [3]; интегрирование уравнения (4) дает 


Г 
Бы 1 (а 60) ЧР ах. (5) 


.14 


Подставив (2) в (5) и исключив из (3) и (5) 
неизвестную величину Г, придем к соотношению 
вида (1) для скорости. | 

Чтобы упростить интегралы (3), (5), восполь- 
зуемся одним общим свойством уравнений со- 
стояния, полученных в [5]: существует давление 
Ре, которое соответствует формированию плотно 
упакованной решетки сфер 


Р, = Раз (о), 


где 0,74 для сфер. 

Если общее давление Р, прикладываемое к 
ДМ, <Рь то даже нижний, наиболее сильно 
сжатый слой мембран не достигает плотной упа- 
ковки (рис. 2, а), и в интегралах (3), (5) удобно 
перейти от интегрирования по координате к 
интегрированию по объемной доле 


(6) 


не аР 
НЙ 9 в 
М=— 0 | не) Чет аа, (7) 
в(Р) 
к ЧР. 
и. =— Со а 4 (8) 
а(Р) 


Первое из этих соотношений можно использовать ' 
для вычисления скорости р 


1 @(Р) ЧР. 
и = Маз ) о 1 (<) те 4, (9) | 
т/п 


второе — для определения толщины мембраны 
<(Р) 
а3‚М | 


Ятп 
"7 ИА 


| ол (©) ЯРаь да 


[ла (а) аа аа. 


(10) 
4% 


Если же давление >Рь, то часть нижних 
слоев мембраны упакована плотно (рис. 2,6) и 
формулы, аналогичные (7—10), примут вид 


м во [60 Р—Р— ; 


ази 
Чит 
РУ \ Г. од га да р (11) 
бо 
бп 
АР. 
ИГ == [а (6%) (Р — Р‚) — | Гл (@) р 4“, (12) 


о 
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И - да [9 —РО+ 
++ \ ОЕ аа, (13) 
тп 
С ЧР 
[ы@9@Ф-—РЭ- | аа] 
Г, = Маз г . 
| Г ЧРаь 
[ль во Р-РЭЧ [ан аа 
лип 
(14) 


Обсудим теперь происхождение и физический 
‘смысл минимальной объемной доли частиц в 
ДМ ашь, фигурирующей в (7—14). В [2] где 
впервые был предложен описанный способ ра- 
счета ДМ, диффузная обкладка мембраны счи- 
талась простирающейся до бесконечности, при 
этом под п (х) понимался избыток числа частиц 
в ДМ по сравнению с их числом в исходном 
растворе. Такая трактовка страдает некоторым 
противоречием, мешающим вычислению именно 
гидродинамической проницаемости: в ней не де- 
лается различия между частицами в растворе 
(которые движутся вместе с ним и не вносят 
вклада в гидродинамическое сопротивление) и 
частицами в мембране (которые не движутся 
вместе с раствором). Понятно, что хотя первые 
частицы и образуют дисперсию с определенной 
объемной долей, но формулы из [3] для гидроди- 
намической проницаемости концентрированной 
дисперсии к ним не относятся. Таким образом, 
вопрос о локализации верхней границы ДМ не 
может быть снят путем отнесения этой границы 
на бесконечность. Толщина мембраны — это ре- 
альная физическая величина, предназначенная 
для деления раствора на две качественно раз- 
личные области: 0<х= Г — ДМ, в которой час- 
тицы не участвуют в нормальном движении жид- 
кости (но могут участвовать в тангенциальном 
[2], не рассматриваемом в данной работе); х>. 
>> — исходный раствор, в котором частицы 
движутся вместе с жидкостью по нормали к 
мембране. 

Величину Г. при помощи формул (10), (14) 
можно связать с параметром &@шш, Который, та- 
ким образом, характеризует состояние диспер- 
сии на верхней границе ДМ; следуя [2], можно 
Положить Яшш равным объемной доле частиц в 
исходном растворе. Можно предложить и другой 
способ определения верхней границы мембраны: 
считать, что на мембране падает определенная 
значительная часть (обозначим ее В) полного 
давления, а оставшаяся часть давления падает в 
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Рис. 3. Зависимость скорости фильтрации от давления при 
постоянном заряде частиц ха==2 (1) и ха=10 (2). 


переходной области между раствором и мембра- 
ной, т. е. там, где частицы движутся с ускоре- 
нием относительно жидкости и тормозятся (в 
этой области, следовательно, не выполняется ус- 
ловие (2) из [2], и нет оснований распространять. 
интегрирование в (3), (5) на нее). Как видно, 
вопрос о локализации верхней границы мембра- 
ны требует специального изучения, поскольку 
предложенные способы только частично отра- 
жают физику процесса. 

В дальнейших вычислениях будем использо- 
вать три способа локализации: 

1) фиксация объемной доли 


а (Г) = би, (15} 
2) фиксация части давления 
а (Г) = &(Р(Г)) = «(ВР), (16) 
3) комбинированный способ . 
(Е) = ши (@шь, (ВР), (17) 


где в8=1—ЗУсшш. Условие (17) при больших 
бита (Опт >10-3) практически эквивалентно 
(16), в котором В=0,9, при малых @шь (17) 
эквивалентно (15). 

Результаты вычислений. Ограничимся изуче- 
нием гидродинамической проницаемости ДМ. 
Подставим в (9), (13) уравнения состояния кон- 
центрированной слабозаряженной — дисперсию 
(соответствующие формулы для случаев посто- 
янного б-потенциала постоянного заряда час- 
тиц приведены в [5]) и вычислим указанные ин- 
тегралы на ЭВМ. 

Результаты расчетов приведены на рис. 3, 4 
(давление и скорость выражены в произволь- 
ных единицах). Отметим прежде всего одну об- 
щую особенность всех приведенных характерис- 
тик: они сублинейны при Р<Роь и линейны при. 
РР. Такое поведение гидродинамической 
проницаемости совершенно естественно следует 
из заложенной в основу расчетов возможности 
сжатия ДМ, ее упругости. При возрастании дав- 
ления ДМ уплотняется и при неизменном коли- 
честве частиц каждая из них вносит возрастаю- 
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Рис. 4. Зависимость скорости фильтрации от давления при 

постоянном &-потенциале частиц и условии обрезания по 

(15) при авиа ==10-8и ха=0,05 (1); 1 (2); 10 (3) (а) и 

условии обрезания по (17) при ха==0,5; аи, ==10-6 (1); 
27.10-8 (2); 10-3 (3) и 8:10-3 (4) (6). 


щий вклад в общую силу трения, поэтому ско- 
рость растет не пропорционально давлению, а 
медленнее. При достаточно больших давлениях 
(когда ашв Уже более не зависит от Р), как 
видно из (13), добавка к скорости прямо про- 
порциональна добавке к давлению, а диффе- 


ренциальная проницаемость Татар 16903" Это 


происходит потому, что при Р>Рь на рыхлом, 
неплотно упакованном слое ДМ падает неболь- 
шая часть полного давления, а плотный слой не 
меняет своей структуры. Таким образом, рас- 
сматриваемая модель позволяет, во-первых, 
предсказать нелинейность гидродинамической 
проницаемости; во-вторых, указать общий ха- 
рактер этой нелинейности; в-третьих, устано- 
вить связь между упругими свойствами концен- 
трированных дисперсий и гидродинамическими 
свойствами ДМ. 

Рассмотрим более детально рис. 3, на кото- 
ром представлены свойства ДМ, в предположе- 
нии постоянства заряда частиц. Интересно от- 
метить, что зависимости на рис. 3 совершенно 
нечувствительны к выбору условий (15—17) и 
очень мало чувствительны кзначению констан- 
ТЫ обши. Это говорит о том, что характер кри- 
вых в этом случае не должен измениться и при 
переходе от условий (15—17) к более реальным 
условиям. С утоньшением двойного слоя (уве- 
личением ха, где х — обратный дебаевский ра- 
диус) давление Ро снижается и графики смеща- 
ются влево по оси Р. Область наиболее ярко 
выраженной нелинейности (заштрихована на 
рис. 3) весьма широка, занимает более двух де- 
сятичных порядков давления (при этом ско- 
рость может измениться всего в 5 раз). й 

Обратимся теперь к случаю постоянного 5-по- 
тенциала частиц. Проницаемость здесь облада- 
ет рядом особенностей: во-первых, даже качест- 
венный ход кривых чувствителен к условиям об- 
резания (15—17); во-вторых, характеристики 
линейны при низких давлениях, что связано с 
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очень резкой зависимостью расклинивающего 
давления от объемной доли [5]. При &=с01$ 
частицы с двойным слоем ведут себя как мяг- 
кие шары и при относительно низких давлениях 
ДМ представляет их плотную упаковку (толщи- 
на которой, как уже отмечалось, пропорцио- 
нальна давлению), в-третьих, область сильной 
нелинейности сравнительно узка (<1,5 деся- 
тичных порядков давления). 

При использовании условия (15) (рис. 4, а) 
общий характер нелинейности лохож на случай 
постоянного заряда, однако с уменьшением тол- 
шины ДС роль нелинейности быстро падает 
(так, при ха >10 расклинивающее давление уже 
практически не сказывается на свойствах ДМ). 

При использовании комбинированного усло- 
вия (17) картина резко меняется: при малых 
значениях ашш характер кривых не меняется в 
соответствии со сказанным о связи условий (15) 
и (17); при больших ошшь появляется падаю- 
щий участок, т. е. область, где Га отрицательна. 
Глубина падения повышается с ростом @шь И 
уменьшением ха; при ха=0,05, ашш=8. 10-3 
(рис. 4,6, 4) и при возрастании давления в 
6 раз скорость падает в 50 раз. Такое поведе- 
ние ДМ связано с конденсацией дисперсии при 
б=с0п${ [5]: при незначительном увеличении 
давления дисперсия резко меняет свою структу- 
ру, превращаясь из плотной упаковки мягких 
сфер (которая соответствует весьма малой объ- 
емной доле частиц) в плотную упаковку твер- 
дых сфер — частиц. При этом гидродинамиче- 
ское сопротивление возрастает очень резко, что 
приводит к падению скорости. | 

Заключительные замечания. Эксперименталь- 
ное исследование обсуждающейся выше квази- 
стационарной гидродинамической проницаемос- 
ти ДМ пока специально не проводилось. Несом- 
ненно, однако, что в свете новых качественных 
результатов, полученных в данной работе, 
представляют значительный интерес для не- 
скольких направлений коллоидной химии и тео- 
ретических основ водной технологии такие ис- 
следования. 

1) исследование структуры 
ния на их функционирование; 

2) исследование релогических свойств кон- 
центрированных дисперсий (в принципе воз- 
можно обобщение результатов данной работы для 
учета вязких свойств дисперсии, сил притяже- 
ния между частицами и, следовательно, пласти- 
ческих ее свойств и т. д.); 

3) исследование строения и режима форми- 
рования двойного слоя частиц, поскольку гидро- 
динамическая проницаемость, как и электроос- 
мос, весьма чувствительны к режиму формиро- 
вания двойного слоя (последнее утверждение 
обосновано в [6]). Поскольку в предельных 
случаях постоянного заряда и постоянного по- 


ДМ и ее влия- 
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Взаимодействие мелких частиц 


верхностного потенциала частиц указанные эф- 
фекты резко различаются, то имеет смысл обоб- 
щение развитых теорий для различных проме- 
жуточных режимов формирования двойного 
_ слоя, например постоянства константы диссо- 
циации ионогенных групп. 

Возможность подобных экспериментов дока- 
зывается результатами, приведенными в [1]. 
Исследовали зависимость скорости фильтрации 
жидкости через гетерогенную ацетатцеллюлоз- 
ную мембрану от давления в области до 2.107 
Па/м?. В результате получались зависимости, 
качественно сходные с приведенными на рис. 
4, а, б (т. е. наблюдалось падение скорости с 
ростом давления). Автор исследования справед- 
ливо связывает такое поведение скорости с ре- 
лаксацией пластических деформаций мембраны, 
которое не имеет прямого отношения к рассмат- 
риваемому здесь механизму, да и сама иссле- 
дуемая система вовсе не похожа на.динамиче- 
скую мембрану. Поэтому ценность цитируемого 
результата для данной работы состоит в том, 
что доказана возможность существования ха- 
рактеристик, подобных приведенным на рис. 
4,6, которые на первый взгляд представляются 
парадоксальными. 

Дополнительные возможности в исследовании 
структуры ДМ представляет сочетание измере- 
ний скорости фильтрации и электропроводности 
ДМ в зависимости от приложенного давления, 
поскольку электропроводность дисперсии весьма 
чувствительна к объемной доле. Некоторую 
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трудность здесь может представить учет поверх- 
ностной проводимости, роль которой возрастает 
при уплотнении дисперсии [7]. Однако несом- 
ненно, что качественные изменения структуры 
ДМ проявятся в ее электропроводности и из- 
мерение последней послужит источником допол- 
нительной информации о динамических мем- 
бранах. 
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ВЗАИМОДЕЙСТВИЕ МЕЛКИХ ЧАСТИЦ С ГАЗОВЫМИ ПУЗЫРЬКАМИ 
ПРИ НАПОРНОЙ ФЛОТАЦИИ 


Т. 3. Сотскова, Ю. Ф. Баженов, Г. А. Голик 


В последние годы напорная флотация, как 
один из методов флотационной водоочистки, 
является объектом постоянного изучения [1— 
3]. Это обусловлено тем, что данный метод в 
применении к очистке воды от высокодисперс- 
ных примесей, образующих устойчивые колло- 
идные системы, обладает рядом преимуществ 
по сравнению с другими флотационными мето- 
дами. К указанным преимуществам относится, 
прежде всего, возможность получения мелких 
пузырьков, размер которых составляет 30— 
120 мкм, стоксовский режим их всплывания и, 
как следствие этого, отсутствие турбулизации 
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жидкости, приводящей к разрушению флотоаг- 
регатов и соответственно к снижению эффектив- 
ности флотационного процесса. 

Большинство работ по исследованию напорной 
флотации посвящено изучению проблем чисто 
технологического характера, сам же механизм 
процесса и влияние физико-химических парамет- 
ров флотационной системы на его эффектив- 
ность изучены в гораздо меньшей степени. Из- 
вестны различные точки зрения, касающиеся ме- 
ханизма напорной флотации [3, 4]. Согласно су- 
ществующим представлениям, формирование 
флотокомплексов может осуществляться за счег. 


и 


